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We study spectral properties of a one-dimensional Dirac equation with various disorder. We use 
replicas to calculate the exact density of state and typical localization length of a Dirac particle 
in several cases. We show that they can be calculated in any type of disorder obeying a Gaussian 
white noise distribution. In particular, we study the random electric potential model, as well as 
a mixed disorder case. We also clarify the supersymmetric alternative derivation, even though it 
proves less efficient than the replica treatment for such thermodynamic quantities. We show that 
the smallest dynamical algebra in the Hamiltonian formalism is u(l, 1), preferably to u(n,n) in the 
replica derivation or u(l, 1|2) in the supersymmetric alternative. Finally, we discuss symmetries in 
the disorder fields and show that there exists a non trivial mapping between the electric potential 
disorder and the magnetic (or mass) disorder. 
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I. INTRODUCTION 

We plan to study some spectral properties of the one-dimensional Dirac Hamiltonian describing a Dirac particle of 
mass m in a magnetic field $ and in an electric (or scalar) potential V: 

h = -ia z d x + <&a x +Ma y + V (1) 

where the fields are random functions of the position x. 

A. A multiple disorder model 

We present a first motivation for this work. In Q the authors drew what the space of coupling constants and 
the renormalization group flow would look like for the 2D random Dirac Fermions. Using bosonization and the 
renormalization group, they derived exact results on localization/delocalization in random magnetic field, and in 
random mass field. However, in the sole presence of an electric potential disorder, the renormalization scheme 
flows away from the Gaussian fixed point to some inaccessible quantum Hall-like transition fixed point. What is so 
interesting with this type of disorder is that magnetic or/and mass disorder do not suffice to account for the integer 
quantum Hall effect fixed point. Thus the electric potential disorder seems to be a relevant field in the flow to the true 
quantum Hall fixed point, even though there are just clues this fixed point exists in the flow diagram of the 2D Dirac 
Fermions. Chalker and Ho in || also emphasized the importance of the electric potential disorder. They took the 
continuum isotropic limit of the network model, known to capture the key features of the localization/delocalization 
transition in the integer quantum Hall effect |lC{] . They managed to map this continuum limit to the model of 2D 
Dirac Fermions in random fields. In particular they showed there exists an electric potential random field locally 
measuring the magnetic flux per plaquette in the network. 

Though it should not shed much light on the 2D problem, we thought it would be interesting to solve thoroughly the 
ID problem. In particular, we shall see that there exists a non trivial exact mapping between the magnetic disorder 
and the electric disorder. This is to be compared to an argument proposed in [g| relating the random mass field and 
the random electric potential field. 

B. Localization 

A second motivation for this work is to know whether a true one-dimensional Dirac particle may be in delocalized 
states. It is already known that there is a zero energy delocalized state in random mass field (or alternatively in 
random magnetic field). Is there any other delocalized states when different types of disorder are present? As we 
shall see, the answer is essentially no. 
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C. From 2D Dirac equation with quasi-unidimensional disorder to the one-dimensional Dirac equation 



The third motivation comes from the study of spectral properties of the quasi-unidimensional 2D Hamiltonian, 
describing a massive Dirac particle in a potential vector and in an electric potential 

h = -ia^ - (T^An - M a z + V. (2) 

The {p, — x,y,z) are the usual Pauli matrices. The potential vector A, the mass M, and the electric potential V 
are all random functions obeying independent Gaussian laws. The three random fields are supposed to depend only 
on x. 

Because of gauge invariance the component A x can be chosen to be zero. Fourier transforming on the y coordinate 
(and thus acting on wavefunctions of the form if>{x, y) = exp(—iky)t()(x)), one gets 

h = -ia x d x - (k + ®)a y - Ma z + V (3) 

where $ stands for A y (x). 

We now fix the average value of M to be m and set k to be the average value of $ thus redefining it into a random 
function with vanishing mean. The average value of V can be chosen to be zero up to rescaling the energy of the 
Hamiltonian. By means of a unitary transformation which applies a x to a z , a y to — a x , and er z to — a y , and for 
convenience, we shall adopt the following equivalent form of the Hamiltonian 

h = -icr z d x + $a x + M (J y + V (4) 

with the given average values ($) = <fi, (M) = m and (V) = 0. The fields $, M and V will be (somehow arbitrarily) 
respectively called magnetic field, mass field and electric potential field, h is precisely the Hamiltonian we would like 
to study. Note that the random mass model (i.e. $ = V = 0) which stems from a particular form of h has already 
been studied in M and H- 



The remainder of this article is organized as follows. In sections P L _0 and III, we map the resolvent calculation 
to a quantum mechanical problem, using the replica trick. In section IVT we derive exact formulas for the resolvent. 
In section [y], we calculate the exact density of states and the typical localization length in various configurations of 



disorder. Then, in section VI, we sketch an alternative derivation using supersymmetry in a straightforward fashion. 



Finally in section VII we study mappings between different disordered Hamiltonians. 



II. MAPPING TO QUANTUM MECHANICS 

We wish to calculate the resolvent of the Hamiltonian h in order to extract its density of states and its localization 
length. In order to do so we will map the associated one-dimensional statistical field theory onto a quantum mechanical 
problem using Feynman procedure. 

A. Path integral representation of the resolvent 

The two-point resolvent is G(x,y,e + if]) — tr 3pin G Q! /3(x, y, e + irf), where a and (3 are spin indices. The spin 
components of the resolvent are 

G a .f3{x,y,e + irj) = (x,a\- \y,0) (5) 

h - (e + irj) 

where r\ is strictly positive. 

Taking advantage of the translational invariance of the ensemble average of the resolvent G{e-\-irj) = (G(x, x, e+ir])), 
one can write 

G(e + ir))= lim ^- [ dx(tr epin G a> 0(x, y, e + irj)) (6) 

which also define G(e + irj) as the average of the previous resolvent G(x, x, e + irj). The integrand can then be formally 
represented by means of a Bosonic field £ with two spinorial components (£t,£x)' so that 
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.. 1 f L , /+ JV^*C a (x)^(x)e X p(i^ L dxdyC(x)(x\h-(e + i V )m(y)). 

G(£ + irj)= lim — — ■ / dx(tr apin r ). (7) 

l^oo 2iL J_ L J ex p (i J_ L dxdy£*(x)(x\h - (e + 

We then obtain 



G(e + if,) = ^ ^ (In I y cxp (t | ^ da** (x) (x| ft - (e + if,) |x)£(x) ) J ) . 



(8) 



The logarithm can be averaged using the replica trick. So replicating n times the two-component bosonic field, we 
are led to 

G(e + i7?)= lim ±Anm|-( fp^rexp^ ^ £ £(aO<s|/i - (e + ir,)\x)t k (x))). (9) 



B. Imaginary time path integral 

Feynman formalism can be applied to map this one-dimensional path integral (statistical field theory) to a 0- 
dimensional quantum field theory, that is to say a mere quantum theory. We shall demonstrate that the corresponding 
physical Hamiltonian is Hcrmitian in some precise sense. Here physical means retaining r\ = in the Hamiltonian 
since the imaginary part of the energy is a mathematical artifice to probe the analytic properties of the resolvent. 

Nevertheless when we re-establish r] ^ in order to obtain analytic properties of the resolvent, we shall rewrite the 
partition function as one of a statistical field theory in imaginary time: 

T>eD?exp{- dxY,£(x){x\-ih+(-ie + ri)\x)S k (x))) (10) 
J - L k=i 

We can now rename — ie + r) as e and make it strictly positive (then shrinking the previous — ie part of the energy 
term to zero). When all calculations are done, we will continue on e to re-establish it to the imaginary axis. From 
now on, G(e) stands for G(e + it]) after having continued — ie + r\ to the imaginary axis. 

As a consequence, we have 

1 P) Pi 

G(e) = -1 lim — — lim ^-Z L {n). (11) 

The i = ^/^T is placed here to compensate the fact that the first derivative is now done with respect to the imaginary 
energy. 



C. Particle- hole transformation 



Actually the previous path integral formulas are ill defined for the reason that they correspond to an unbounded 
quantum theory. To see it, we follow the Feynman procedure, write the classical Hamiltonian associated to the 
previous path-integrals, then quantize it. The momenta for the spinorial fields are 

n> = £ = -fa* (12) 

where C is the Lagrangian density and the dot indicates derivative with respect to x. So the Euclidean coherent state 
path integral formalism yields the classical Hamiltonian 



n ~ 
7, — 1 J 



(13) 



fe=i 1 

before any ensemble averaging. To quantize this Hamiltonian, we impose the following canonical commutators: 

My),U il (x)}=6 k , l 6(x-y). (14) 
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This leads to a e.H^ k a z ^k term in the quantum Hamiltonian which makes it unbounded. This is due to the minus sign 
appearing within the spin down components product. To quantize the Hamiltonian safely wc have the freedom to make 
any transformation on the canonical operators whenever the canonical commutation relations remain unmodified. For 
Dirac Fermions the correct procedure would be to fill the Dirac sea upon making a particle hole transformation. 
However, such a change of vacuum is not appropriate for the Bosons that would condensate anyway. So healing the 
unbounded Bosonic Hamiltonian requires a stronger transformation. One can set 

n£ fc = (&f,-&i) & = (& T ,&f)- (15) 

The minus sign is due to heal the unboundedness of the spectrum. The additional particlc-holc transformation 
restores the commutation relations destroyed by the inversion of the spin component operators. The tribute to such 
a regularizing transformation is that the spinor momentum n^ fc is not any more the Hcrmitian conjugate of We 
shall see that the internal dynamical symmetry of the replicated Hamiltonian is not U (2n) any more but the non 
compact U(n, n) group. We will call dynamical algebra (or improperly group), any algebra which generators are the 
sole building blocks of the Hamiltonian. 

Let us define the Hermitian spin operators Jki = (icr x , icr y , <j z )£,i- With the operators Qm = in^ fc ^, they form 
the 4n 2 generators of u(n, n). In components, the spin operators read: 

J& = £(&W-&T&l) J « = ^(&t&J"+M'l) J« = ^(&f&T+&t&i)- (16) 

We will notice that, be it before of after averaging, the Hamiltonian (or the Hamiltonian density) does only depend 
on the sum operators 

n n n 

fc=i fc=i fe=i 

Then we define the corresponding ladder operators, which we shall later use to generate physical sectors of the effective 
Hamiltonian 

Z+ = Z V - iZ x Z~ = Z y + iZ x Z° = Z z . (18) 

They obey su(l, 1) commutation rules 

[Z+, Z~] = -2Z° [Z°, Z~\ = Z+ [Z°, Z~\ = -Z- (19) 

which have to be completed with the identity (Z + Y = Z~ , which tells us which real form of sl(2, C) the spin operators 
Z are the generators of. 

This statement implies that the dynamical algebra can be reduced from u(n,n) to u(l, 1). Such a strong simplifi- 
cation, which was not previously known, also occur in the supersymmctric version of the derivation. 



D. Averaging 

The disorder fields can be written away in the vector field A = (M, <£, V) with mean value A = (m, 4>, 0). The most 
general Gaussian disorder distribution can be written using a covariance matrix C. C is assumed to be symmetric 
and positive definite, as it seems physically reasonable. The Gaussian measure is then 



(20) 



P(A) = cxp (-1 y°° (A - A„)t<7(A - A )rfA^) . 

Performing the average, one obtains an effective Hamiltonian of the form 

iff = -A .Z - Z^C^Z (21) 

which does not depend any more on x. We next rescale our Hamiltonian in order to eliminate the inconvenient | 
factor. Choosing to discard off-diagonal mixing terms in the covariance matrix, we will adopt the definitive form of 
our effective Hamiltonian: 

H = eZ° - imZ x + + g x {Z x f + g^f + g z {Z z f (22) 

where g x , g y and g z denote the disorder strengths. As announced, the effective Hamiltonian does only depend on 
the spin operators Z. From the definition we retained for a dynamical algebra, one immediately sees that u(l, 1) is a 
dynamical algebra of H. 
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III. GROUNDSTATE SECTOR 



We shall now describe some groundstate sector properties in a view to evaluate Z^ (n) when L goes to infinity and 
n goes to 0. 



A. Basis description 

One can look for the groundstate of the effective (averaged) Hamiltonian in the lowest weight space generated by 
the dynamical symmetry w(l, 1). Since the trace operator of u(l, 1) commutes with H, we can restrict our dynamical 
symmetry to su(l, 1) in any labeled sector of the Hamiltonian. The lowest weight of the desired space is equal to n/2 
as the sum of the n spin-1/2 discrete lowest weight representations of each individual spin. Call |f2) the lowest vector 
of this particular representation. The states of this representation space, generated by the action of operator Z + on 

|fi) are \k) = » k / |0). The norms of such states can be quite easily computed. We obtain 

Reducing the dynamical symmetry allows us to generate the whole representation space as above. Another advantage 
is that one can as easily calculate the action of the spin operators in the {\k}} basis. 

Z°\k) = (k + Z + \k) = (k + l)\k + l) Z~\k) = (k + n-l)\k-l). (24) 



B. Interpretation in terms of representation 

We shall see now that the replica limit n goes to is mathematically well defined. As long as n is a positive integer, 
the groundstate sector remains a lowest weight representation space of su(l, 1) belonging to the discrete series. When 
n tends to in the interval ]0, 1[, then the representation space does not any more belong to the discrete series but to 
the supplementary series. In Vilenkin's notation, this representation is labeled by (I = n — 1, e = 0) The series 
of \k) states is not any more truncated below by the action of Z~ on |1 — n). Theoretically it spans down to — oo. 
Remarkably, when n goes to 0, we have 

d*>-f^5TiJ^^- (25) 

So that, at the precise point corresponding to the physical disordered system (n = 0), the restriction of the scalar 
product within the groundstate sector degenerates. The \k) states are null states except when fc = 0. When n is in 
]0, oo [, one can define a proper scalar product in the groundstate sector for wavefunctions of the form \<p) = ^2 k=0 4>k\k): 



which is identical to Vilenkin's formula (see also juj) but truncated below k = 0. 



IV. GROUNDSTATE ENERGY OF THE REPLICATED HAMILTONIAN 



A. Groundstate wavefunction equation 

Recall the effective replicated Hamiltonian is 

H = eZ a - imZ x + icj)Z y + g x (Z x ) 2 + g y {Z y f + g z {Z z f . (27) 

The expansion of the groundstate on the {\k)} basis is |^) = X^fclo Ck\k). It satisfies the eigenvalue equation H\^f) = 
A(n)l'J'). A(n) is expected to behave like An + 0(n 2 ) when n goes to 0. Re-written in {\k)} the basis, this eigenvalue 
equation yields 
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«k + -)+ g z (k + -) 2 )C fe - -(fcCfe-i - (* + n)C*+i) + *|(*Cfc-i + (* + «)Cfc+i) 

1 1 n 1 

+ -fc(fc - - ff x )Cfc-2 + -x(k(k + n) + -){g y + g x )( k + ^k(k + l)(g y - g x )( k +2 = A(n)Ot 

The equation is valid for any integer k > 0. 

B. Left and right eigenstates of the Hamiltonian 

For technical reasons, that is when rj ^ 0, we must consider a non-Hcrmitian Hamiltonian. As a consequence, 
Feynman's formula which gives probability amplitudes of a quantum mechanical theory in terms of one-dimensional 
statistical path integrals, must be written in terms of left and right eigenstates. This follows from the resolution of 
unity used to derive the Feynman equivalence. Indeed, the resolution of the identity is expressed in terms of right 
and left eigenstates. 

A simplifying feature is that our Hamiltonian and its Hermitian conjugate are unitarily related by means of the 
simple unitary transformation R — exp(i7r X)fc=i ^kk) = exp(i7rZ 2! ). Left eigenstates lie in the dual vector space of 
the current Fock space. In principle, there is no need to project them to the original vector space. But doing so is 
convenient. Using the transformation R, we can map those left eigenstates back into the Fock space: 

\lp)l = R\1>)r. (28) 

C. Groundstate wavefunction 

In order to apply Feynman formalism, we must normalize our groundstate wavefunction. When n ^ 0, we have 

In the limit n goes to zero, we obtain || , 3 / || 2 = |Co| 2 - So that the normalization condition can be chosen (o(n = 0) = 1. 

D. Hierarchical equations 

Since we are interested in the limit n goes to 0, we expand the groundstate wavefunction and its eigenvalue in 
powers of n (see Ol for the use of such scheme on the Halperin model). At order 1, we set the notations 



A(n) = n\ + 0(n 2 ) 
Ck(n) = ip k + nxk + 0(n 2 ) 



2^ (30) 



At order 0, we find the equation for the groundstate of the physical system (i.e. when n = 0) 

(e + {g + 9x * 9v )k)i/j k ~ y (V'fc-i - 1>k+i) + ^(V'fe-i - i>k+i) 

+^(9x - g v )(k - i)ipk-2 + ^(g x - g v )(k + 1)^+2 = 0. (31) 

The result is an extension of an equation obtained by Luck [Q , then Balents and Fisher [^) . The coefficients ipk should 
bear a physical interpretation in terms of amplitudes of multiple scattering on impurities. At least formally, in the 
random mass model with zero mass average, the coefficient ip2k may be interpreted as the right-hand-side Berezinskii 
block R k (see [@, 0). The recursion relation on the Berezinskii blocks R k or on the groundstate coefficients ip2 k 
are indeed the same (the odd V^n+i do not contribute). Similarly the coefficients of the left groundstate may be 
interpreted as the left-hand-side Berezinskii block R k ■ Those blocks are the sum of all contributions due to scattering 
events happening on the right and on the left of the point where the one-point Green function is computed. 

Coming back to the hierarchical system, order 1 yields the equation for Xk'- 



G 



(e + (.9 H ^ )k)Xk - -j(Xk-i - Xk+i) + i 77 Ufc-i - Xfc+i) 

+ |(.9x - - l)Xfc-2 + ^(ffz - 9y){k + l)Xk+2 

+ + ~)^fc + gktpk + ^( 2fc + - 9x)^k+2 + ^(m + i(f>)ip k +i = A^fe- 

The equation at order 1 can then be legitimately evaluated at A: = and yields the surprisingly simple result: 



9x + Qy . m + i<f>ipi g y - g x fa 



(32) 



E. Resolvent formula 

From formula (jl^) it is sufficient to evaluate the partition function Zl (n) in the limit where L goes to infinity and 
the replica index n goes to zero to calculate the imaginary time resolvent. In that limit, we have 

Z L (n) ~ exp(-if)|*) r - exp(-2LAn + 0(n 2 )), (33) 

so that with equation ( |ll| ) 

G(e) = 2z — . (34) 
Note that the 2 factor is due to the adopted scaling convention of the Hamiltonian, and consequently of A. 

F. Sum formula 

In the previous subsection, we derived an efficient way to obtain the resolvent. Now, we shall show that though it is 
more cumbersome, we can also write the resolvent as a sum as was done in At first, let us replicate the correlation 
function n times without changing its value: 

1 JjZWj ELi ex P (- JZo dx E n k=i gfeKj ~ih+ He + r,)\x)M*)) 

(€) i (t JVm*^p(~J^d X j: n k=i e k (x)(x\- l h+(- l e + r 1 )\ X ^ k (x)) h [ > 

In the limit where n goes to 0, the partition function denominator goes to 1, so that averaging becomes possible 

1 p 1 71 p CO n 

G(e) = - lim (tr spin / V&>?- £ £*{x)$(x) exp (- / dx^Ckix^l -ih+ (-ie + r,)\x)^ k (x))). (36) 

After regularizing this correlator by particle-hole transformations as was done above, and mapping the path integral 
to the equivalent quantum problem, we get 

G(e) = -2ilim-(V\ l Z m r . (37) 

To render explicit this amplitude, we make use of our scalar product 

G(e) = -2i lim ^f)(-l)»(fe+ g)l^i 2 Jtt^Y (38) 
Eventually taking the limit n — > we obtain: 

oo 

G(e) = -2iJ2(-V k m 2 . (39) 

k=Q 
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V. DENSITIES OF STATES 



We now contemplate various cases of disorder and calculate the subsequent resolvents. 



A. Random mass model 



The random mass model, which corresponds to the original Hamiltonian h = —ia z d x + Ma y is the only previously 
(well-) known case. M is the random mass field and follows the Gaussian measure P(M) = exp(— ^ dx(M-m) 2 ). 
It is equivalent to many one-dimensional physical models (see 0j for a review on the quantum supersymmetric 
Hamiltonian ft 2 ). 

After averaging the effective Hamiltonian one obtains 



H = eZ° — imZ x 
Correspondingly, its groundstate equation reads for k > 1 



(e+ -k)ip k - y(V>fc-i 



whereas the groundstate energy (32) is slightly reduced to 

A = 



9 
4 



m tpi 



4</V 



(40) 



(41) 



(42) 



Such a difference equation can be solved exactly. In literature the technique of the generating functional is often 
used to solve, when possible, quasi-linear recurrence equations. We prefer to use Laplace method which, though 
equivalent, is much more convenient in such circumstances. Indeed it allows to pick up a particular solution of the 
equation thanks to the choice of the appropriate contour in the complex plane. This is all the more useful since we 
expect only one physically acceptable solution among a finite dimensional vector space of solutions. 

Let us detail how we use it. We start giving an integral representation of ipk in the complex plane along a contour 

C: 

i\>k = J dwK(k,w)ip(w). (43) 

The kernel K has to be chosen in order to yield the simplest differential equation on il>{w). The contour C must 
be chosen so as to obtain a solution different from zero. It also has to be chosen in a view to pick up the desired 
solution in case there are more than one. If the contour possesses boundaries, we require ip w to tend to zero fast 
enough (depending on K) in their vicinities. 



For our purpose, we choose the Melin kernel, K(k,w) 
correspondence rules 



,-fc-i 



Then it is straightforward to show that we get the 



kif^k — * wd w tp u 



Tpk+l > 



w 



W1p u 



The groundstate wavefunction equation then reads 



m 1 
~2^w 



5/2 1 > 



4>v 



g 



l 



w 







4 

which is a solvable (first order) differential equation. An acceptable contour for which the series 



oo 

E 

fe=0 



(-1)1^ 



(44) 



(45) 



(46) 



(on which depend G(e)) converges is the interval ] — 1, 1[. Solving for ip w then inserting it in the integral representation 
(fl3|), we finally obtain 



ipk 



dw — 
l 1 



exp(- 



g 



-)■ 



(47) 
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It is noteworthy that this integral has already been obtained in Q. In this paper, the author treats the problem of the 
XX spin chain, with disorder in the couplings, our mass disorder, and bathed in a transverse magnetic field, playing 
the role of our mean value m of the random mass field. His derivation resorts to the characteristic exponent f2(e) of 
the localization problem, which can be formally identified to 2A(ie). This model is known to be ultimately equivalent 
to the random mass model with zero or finite mean value of the random mass field. 

The latter ipk is not normalized. Recall \\^\\ 2 = IV'ol 2 = 1, so taking tpk/ipo as the wavefunction component 



normalizes the groundstate. Thanks to the change of variable w = tanh(ij) 
finally obtain an expression for A 



and to some integrations by parts we 



m 

T 



e *i-?(f) 
2Jf_-(f) 



(48) 



where K is the modified Bessel function of the second type. Hence, the resolvent may be put in the form: 



G(e) = -i- 
oe 



Ki 



■(f) 



(49) 



which is the result obtained in |jj f° r h 2 . The prime denotes derivative with respect to the argument of the Bessel 
function. From this imaginary energy formula, we can derive the integrated density of states as well of the inverse of 
the typical localization length. Indeed the application of Thouless formula to our case yields 



2X(ie) = r l (e) +mN{e) 



(50) 



after analytic continuation on the energy. The final expressions for the density of states and inverse localization 
length can be found in §. They are drawn in figure [l]. They are special cases of the foregoing formulas related to the 
problem of a random mass particle in a constant magnetic field (see table ||) . 



B. Random electric potential model 



1. Massless Dirac particle 



The random electric potential model of a massless particle corresponds to the original Hamiltonian 

h = -ia z d x + V (51) 

where the electric potential V follows the Gaussian measure P(V) = exp(-^ dxV 2 ). The model is trivially soluble 
because the wavefunction eigenstates can be explicitly written at all energies 

ip e (x) = exp(ia z [ [e - V(u)]du) (52) 
Jo 

and for any quenched disorder V. 



2. Massive particle 



The random electric potential model of a massive particle corresponds to the original Hamiltonian 

h = —i<J z d x + m<7 y + V (53) 
where m is the fixed mass of the particle and V, the electric potential follows the Gaussian measure P(V) — 

ex p( h f-°o dxy2 ) ■ The effect 

ive Hamiltonian is H = eZ z — imZ x + g(Z z ) 2 . Correspondingly, the groundstate 

equation is for k > 1 

Tfl 

(e + gk)i/j k + yCV'fc+i - V'fc-i) = (54) 
whereas the groundstate energy is reduced to 



9 



e mipi 

A = - H — . 

2 2 



(55) 



The careful reader might have noticed that this recurrence relation is the one for modified Bessel functions. Since we 
must cautiously take into account boundary conditions, it is nevertheless safer to use again Laplace method. For our 
purpose, we choose once more the Melin kernel, K(k,w) — to —1 . The groundstate energy then reads 



to. 1 . 
e+ -(--w) 
2 w 



ip w + gwd w ip w = 0. 



(56) 



An acceptable contour is precisely one of those used for the Bessel functions integral representations. The contour 
comes from — oo slightly below the real axis, rounds the origin then goes back to — oo slightly above the real axis. 
Solving for if> w then inserting it in the integral representation, we obtain the non-normalized value of ipk 

exp( — (w + —)). (57) 
oo 2g w 



We finally get an exact expression for A 



A 



Ah 



■(-) 

v a ' 



9 v a ' 



where / is the modified Bessel function of the first type. Thus, after a little algebra, 

d 



G(e) 



1 a 



a y 



s v a ' 



(58) 



(59) 



In the disorder free model, the mass band ] — m,m[ is a forbidden energy region. The density of states and 
localization length are well-known and given in table 0. 

Coming back to the disordered case, the density of states and inverse localization length can be obtained through 
Thouless formula. In that case, it is cumbersome to write definitive expressions for them because the Bessel functions 
I bear an imaginary index. We have drawn them anyway in figure ^|. 

When there is disorder, the particle penetrates (statistically speaking) in the mass band. That is why it is interesting 
to study the weak disorder limit in that range of energy. To derive asymptotic expansion of the density of states and 
inverse localization length, we apply the steepest descent method to the integral representation of the modified Bessel 
function /. Because the index of the Bessel function is an imaginary index, we cannot straightforwardly resort to 
tabulated asymptotics Bessel function expansions. Surprisingly, this technique only yields real values when expanding 
X(ie). This remark is valid to all orders. The imaginary part of X(ie) that accounts for the part of the density curve 
in the band ] — to, to[ is provided by the direction of rapid oscillations which is transverse to the steepest way in the 
complex plane. That is why Im(A(ie)) cannot be accessed through the steepest descent method. 

The steepest way expansion is worked out within the |e| < to band. As for the inverse localization length, we find 
in the limit g goes to 



d f 1 

l~ l {e) ~ -gm— < - In 
dm [ 2 



,H — (cot(a) + a) > where a = arcsin( — ). (60) 
27TTOCOs(a) / g J to 

As for the integrated density of states, when the energy is smaller enough than the mass, the imaginary part of the 
Bessel function becomes crucial. Fortunately, the imaginary part of the modified Bessel function of imaginary index 
can be accessed directly 



777 

Im(/«(-)) 
9 9 



sinh(- 



2i sin(wr-) 

v a ' 



9 9 9 9 



sinh(-) m 

^s(-) 

7T f g 



which allows an approximate evaluation. We then find 



N(e) 







-gm— — arctan < sinh( — )\Jcos{a) exp( 

dm I g g 2gm g 



(cot(a) + a)) 



On average, the residual densities of states at e = is 



P(0) 



2m 
ng 



(61) 



(62) 



(63) 



The average length of energy penetration in the mass band can be extracted from ( [62] ) and is of the same order 
than p(0). 
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C. Multiple disorder model 



We would like to know what the effect of combining two or three types of disorder is. We will show that such 
a combination destroys the e = singularity. To be as simple as possible, we will choose the mean values of those 
random fields to be vanishing. Yet the calculations could be generalized to non-zero mean values. This corresponds 
to the original Hamiltonian 

h = -ia z d x + $<r x + Ma y + V (64) 

where M is the random mass of the particle, $ the random magnetic field and V the electric potential. 

The effective Hamiltonian is H = eZ°+g x (Z x ) 2 +g y (Z y ) 2 +g z (Z z ) 2 . Correspondingly the groundstate wavefunction 
equation reads for k > 1 

(e + (fh + 9 -^ ]L )k)^k + \{g y - 9 x )(k - l)V*-2 + \(Sv - 9x)(k + 1)^+2 = (65) 
whereas the groundstate energy is 

X= i + 9_x + 9y 9_v^9x^ (66) 
2 4 2 

Only the even terms contribute. Indeed one can show that the odd terms of the sequence {ipk} yield an unacceptable 
diverging sub-series. 

So let us set ip2k — * ipki define S = g y — g x and a = 2g z + g y + g x then the previous formulas turn to 

A _ e _l 9x + 9 y ^ 9 V - 9x ipi 



2 4 2 Vo 

S 1 S 1 

(e + CT.fc)Vfe + -(* - -)4>k-i + 7;( fc + nWh+i = o. 



2 V 2' 



Now introduce the roots of w 2 + 2jw + I = 0: 



and the constant 



a± = --5 ± tp- 1 (67) 

/3=T7— r. (68) 

o(a+ — a_J 

Since | j| = | 2g ^ +9x ^ +gy | > 1, those roots are always real. This is what will guarantee that a critical point at e = 
cannot exist in presence of at least two types of disorder. To make what follows clearer, we choose 6 < so that the 
roots a± are positive. 

Using the Laplace method, solving for ip w , then we find 

ipk = — <b dw.w~ k ~^ .(w - a+)~^~^ e (w - a_)~5 +/3e 

2*^ J[a+,oo[ 

cos(7r/3e) f 00 , (a-a + )-*-' 3£ (a-a_)-i+ /3e 
' da. 



so that we obtain A and 



leJ V + 2 &/~da.a-i.(a-a+)-4-' Je (a-a_)-i+/'^- ^ 



Coming back to the original notations, we have: 



2.9^ + 9x + 9y , If 1g z + g x +g y ^ 2 



a± =- ' a * ' atf ± W JZ ■ ^ ■ ^ _i and fl= — (70) 
5y - 9x V V 5* ~ 5* / 2^^ + 5xK^ + 9y) 

The typical inverse localization length associated to this imaginary energy resolvent is finite. The presence of a second 
type of disorder (or more) makes the random mass critical point massive. The density of states and inverse localization 
length are drawn in figure [| 
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D. Random mass particle in a constant magnetic field 



This model is equivalent to a Dirac particle with a fixed finite mass in a random magnetic field. It corresponds to 
the original Hamiltonian h = —ia z d x + (j>a x + Ma y . (f> is the homogeneous magnetic-like field and M is the random 
mass field which follows the Gaussian measure P(M) = exp(— ^ dx{M — to) 2 ). (f> can also be interpreted as a 
uniform staggered magnetic field in the third direction perturbating an XX random spin chain. After averaging the 
effective Hamiltonian one obtains 



H = eZ" + {(/)- im)Z x + g(Z x ) 2 



(71) 



The calculations are of mino r dif ficulty in comparison to the random mass case provided one uses the non-trivial 
mapping described in section VII. So we just give the results: 



G(e) 



de 



\Je 2 + (jr 



KM 



(72) 



Continuing on e we obtain the expressions for the density of states and inverse localization length given in table 0. 
They are also drawn in figure ^. 

In this section, we have derived imaginary energy resolvent formulas which bring about the exact density of states 
as well as the exact typical localization length in the cases of a few interesting models. As long as the disorders are 
Gaussian white noise, we can calculate those quantities exactly. This is due to the fact that equation ( |3l| ) is linear 
in k. The Laplace method then provides us with a first order differential equation. That is why an exact explicit 
expression can always be obtained. 



VI. SUPERSYMMETRIC ANALYSIS 

In S, the authors make use of the supersymmetric trick to describe the random mass model in the vicinity of 
e = 7Ti = 0. They obtained equation (|4l| ) after they identified the groundstate sector by pre-diagonalizing the effective 
Hamiltonian. We will see that this diagonalization is not necessary and that the sought sector can be algebraically 
generated as a degenerate su(l, 1) grade-star representation. Similarly to the replica version, we will identify U(l, 1) 
as the dynamical group of the effective Hamiltonian, preferably to £7(1, 1|2). 

The path integral representation of the resolvent can be obtained as in g . We recall the formula for the averaged 
resolvent: 



G(e + irf) = (tr.^ / V^VCThpT>i/}.^ a (x)ilip{x) exp(- / dx?(x)(x\ -ih+ (-ie + r))\x)£(x) 

J J — oo 

dx7p(x)(x\ — ih + {—ie + 7i)\x)ip{x))) . 



As in the replica case, we can now rename — ie + 77 as e and make it strictly positive. When all calculations are 
done, we will continue on e to re-establish it to the imaginary axis. Similarly to the replica case, the previous path 
integral formula is ill defined, at least as far as the Bosons are concerned. To see it, we follow the Feynman procedure 
and write the classical Hamiltonian associated to the previous path-integral. We will then quantize it. The momenta 
for the spinorial fields £ and ip are 

IL/i = = -ijja z and 7T£ = ^ = -£*<7g (73) 
dip ot, 

where the dots indicate derivative with respect to x. Here the Fermionic derivative is a right derivative. The Euclidean 
coherent state path integral formalism yields the classical Hamiltonian 

H c = J dxU^.d x ij + C(U i „7p) + U i .d x ^ + C(U i ^) (74) 

before any ensemble averaging. 
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To quantize this Hamiltonian, we impose the following canonical commutator: 

Wy),n^x)} = S(x-y) [Z(y),n 6 (x)] = 6(x-y). (75) 

This leads to a eJI^er^ term which makes the quantum Hamiltonian unbounded, due to the minus sign appearing in 
the product of the spin down component operators. To quantize the Hamiltonian safely we have the freedom to make 
any canonical transformation on the spin component operators. For the Dirac Fermions the correct procedure is to 
fill the Dirac sea upon making the particle hole transformation: 

n^ = (^,Vi) ^ = (^T^ i + )- (76) 

However, as we already noticed in the replica case, such a procedure is inappropriate for the Bosons. As before we 
need a stronger transformation. We then set 

n 5 = (£ T + ,-£i) e = (^T^i + )- (77) 

The minus sign heals the unboundedness of the spectrum and the additional particle-hole transformation restores 
the commutation relation destroyed by the inversion of the spin down component operators. Again, as a result, the 
spinor momentum is not the Hermitian conjugate of £. The internal dynamical symmetry of the supersymmetric 
Hamiltonian is not any more U(2\2) but E/(2|l, 1). 

In order to go beyond this dynamical symmetry of the Hamiltonian, let us define the following Hermitian spin 
operators 

j x = - m) + - wo j v = W^t + mi) + IftW + wo 

Like in the replica case, we point out that the Hamiltonian does only depend on those sum operators. Then we define 
the corresponding ladder operators 

J+ = J y - ij x J- = J y + ij x J° = J z (78) 

which obey su(l, 1) commutation rules: 

[J + ,J-] = -2J° [J°,J-]=J + [J°,J-] = -J-. (79) 

However, as such, this statement is incorrect because we cannot avoid {J + )^ ^ 3~ ■ By computing the commutation 
relations of the spin operators, we merely proved they generate a sl(2,C) algebra but did not stipulate its real form. 
The reason of that difference is that although the Bosonic part of J + represents a Bosonic su(l, 1) spin, the Fermionic 
part represents a su(2) spin. 

Now we shall show how to "transmute" this last spin into a su(l, 1) one. Such a transformation is conceivable 
because the Fock space is a grade space. Hence two types of star representations are permitted. We shall resort to 
the unusual one: the grade-star representation. It is characterized by the properties of the adjoint operation on 
the Lie algebra operators: 

(aP + bQ)^ -SPt+bQt [P,Q]^ (-l) dcsP - dcsQ [P\Q^] (Q+)+ = (-l) dcgQ Q. (80) 

a, b are c-number and P, Q are operators. degP takes the value or 1 whether P is even or odd. 

In order to satisfy such identities, we must change our Fock space scalar product. It was previously implicitly 
defined by the creation/annihilation operators and how they act on the vacuum. Now assume that instead of the 
usual definition of conjugation for Fermionic operators 

(V« T ) t = ^ W»f) f = ^ and {^=4>t (^=^1 (81) 

we prefer 

(V, T )t=^+ (^+)t = _^ T an d tytf = il>t (V+) t = -Vi (82) 

which are coherent assumptions provided we accept to redefine the scalar product in the Fermionic sector of the Fock 
space. One can check the compatibility of this new scalar product with the grade-star adjoint operation through 

(P^x\y) = (-l) d ^ p - dc ^(x\Py) (83) 
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for any homogeneous operator P and homogeneous vectors x and y. Yet, this scalar product is not any more positive 
definite (that is why in principle we should not call it any more scalar product). For example if |f2) is the Fermionic 
vacuum of norm unity, then | \ip* \Q)\\ 2 = 1 and 1 1 2 = — 1- More generally a state with An or An + 1 Fermions 

has norm 1 , whereas a state with An + 2 or An + 3 Fermions has norm — 1 . 

A consequence of this change is that we recover {J + )^ = J~ ■ Another related consequence is that the effective 
Hamiltonian in real time and when r\ = (physical limit) is a true Hermitian operator, provided we accept the 
degenerate scalar product. As a result the J spin operators now thoroughly obey the constraints of a su(l, 1) algebra. 
This statement implies that the dynamical algebra can be reduced from u(l, 1|2) to u(l, 1). Such a strong simplification 
is analogous to what occurs in the replica derivation. 

A practical consequence is that once we identified the (grade) star representation of the dynamical symmetry of the 
Hamiltonian, we can easily generate its main sectors. For example the groundstate sector is generated by the states 
|fc) = where now |fi) stands for both the Bosonic and Fermionic vacuum. Within our formalism we can 

now write them explicitly in terms of creation operators: 

|fc> = = m W+<W)*i°> = h.^^ )k + Me T % + ) fc -VM + ) fc N) (84) 

which are precisely the states obtained in |§. 

The left eigenstates can be mapped from the vector dual space to the usual space by means of the unitary operator 
R — exp(i7r^7 2 ) [], which is a supersymmetric operator. The action of the general Hamiltonian 

H = eJ°- imj x + u/>jv + g x {J x f + g y {J y f + g z {J z ) 2 (85) 

on the expanded groundstate = X^oCfcl^) produces 0-order equation (^l|). Contrary to the replica analog, the 
supersymmetric method only leads to the sum formula ( |39| ) for the resolvent. Indeed, using the mapping from path 
integral formulation to quantum mechanics, we may only write the resolvent as 

G(e) = -2*(tf|,5°|#) r . (86) 



where 5° is the fermionic part of the generator J° . It is the analog of (p7U and leads to formula (39) 



VII. SYMMETRIES OF DISORDERS 

In this section, we discuss some mappings that may exist between Hamiltonians with different disorders. Such 
question is of most interest in the 2D disordered Dirac Fermions model. Indeed it has been argued that the renormal- 
ization group flow of the random mass model can be, in some sense, mapped to the random electric potential model 
P|. In this one-dimensional case, not only does it exist at least one non-trivial mapping between disorders, but it 
can also be made explicit. Besides there exists a more obvious mapping between magnetic and mass disorder which 
deserves quoting. 

The geometric group with double covering group SU (1, 1) is the Lorentz group 50(2, 1). So we may imagine the 3D 
space of disorder coupling strengths in the following way. The time (Oz-)axis measures the potential electric coupling 
strength whereas the Ox and Oy axis measure the mass field and magnetic field strengths. The 50(2, 1)— invariant 
manifolds are the hyperbolo'ids of axis Oz. J x , J y and J z are the generating operators of the Lie algebra. Precisely 
J z generates the Euclidean rotation around axis Oz, whereas J x and J v generate boosts of axis Ox and Oy. As a 
consequence, we can write the action of a half-turn rotation around the time axis on the spin operators: 

e -™j z jx e ™j z = jy e -i*S* jv e ™j° = _J* ( 87 ) 

This rotation operator establishes a unitary mapping between physics of the random magnetic field and physics of 
the random mass field. 

More interesting is the remark that no unitary transformation can map a mass disorder (or a magnetic disorder) 
to a potential electric disorder. Indeed no unitary boost is able to map a spatial axis to the time axis. Let us now 
consider the non-unitary, yet non-singular, operator R = exp(-|J rj ') which is a boost of axis Oy and imaginary angle 



1 It must not be confused with the Bosonic operator used in 
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We will show that it maps a random mass problem to an electric potential model. In particular we have these 
actions on the spin operators: 

e^ JV J z e-^ JV = ij x ei jV J x e~i jV = ij z . (88) 

Let us start reasoning from the effective Hamiltonian H = ej z — imj x + g(J x ) 2 of the random mass model. Then 
we make R act on H : 

RHR- 1 = tej x +mj z - g{J z f . (89) 

Upon making the replacements 

e — ► m m — ► — e and g — > — g (90) 

we turn the random electric potential problem of section [v| into the model describe by RHR^ 1 . 

However being non-unitary the transformation R is not trivial. Above all we cannot use scalar-product-dependcnt 
properties of the Hamiltonians since we lost unitarity. That is why in such circumstances, we must consider our Fock 
space not as an Hilbert space but as a mere vector space. 

Yet we notice that the kernels of operators are covariant objects, i.e. Kei(RHR^ 1 ) = R(Ker(H)) because R is non- 
singular. Nevertheless under suppression of the scalar product, the kernel of H has enlarged from a one-dimensional 
subspace to a two-dimensional one, because any normalization condition has been dropped. Let us first describe 
Ker(e l / 2 — imj x + g{J z ) 2 )- We already know that the wavefunction 7fc_|_|(^) belongs to the kernel. A second 

independent solution is (— l) k Kk+± (y)< To check this last statement, change variable in the recurrence equation (|54|) 

from tpk to (— l) k ipk, and observe it became the recurrence relation for the modified Bessel function of second type. 
So we conclude: 

KeriRHR- 1 ) = Vcc{/ fe _ f {-),{-l) k K k ^ (-)} (91) 

9 g 9 g 

Our task is then to find out in a two-dimensional kernel what the physical image of our original groundstate is. We 
must select the right combination in the 2-dimensional space which stands for the groundstate wavefunction of our 
primary problem to perform such a selection. This combination does not depend on the values of the parameters of 
the Hamiltonian because R acts linearly on the states. 

To proceed we can use the alternative method of the replica and view n as a deformation parameter. Then the 
groundstate energy of the replicated Hamiltonian X(n) = \n + 0(n 2 ) is characterized by 

2 2 

In the supersymmetric version of the problem, the selection rule of the groundstate is the convergence of the series 
defining the correlation functions. It applies in particular to the resolvent. In the replica case, the analogous selection 
rule is the sign of A which states whether the quantum theory is correctly defined (bounded or not). Since we had to 
forget about the scalar product, it is cumbersome to evaluate the series defining the resolvent in the original problem, 
whereas it is much simpler to check the sign of A for any reasonable values of parameters e, m, and g. 

So we can now apply our selection rule. Essentially because a minus sign appears in formula fl9^) in comparison to 
the original formula (pq), the only acceptable combination in Kei(RHR~ 1 ) is the wavefunction ip^ = {—l) k K^m{^). 
Finally, we conclude 



G (e) = 2i|-A = 

oe ae 



f — 9 —^- 

9 y g' 



(93) 



and recover result ([19|). Practically, this way to obtain (|49j) saves many calculations. In particular we transform a 
recursion relation of order 5 to another one of order 3, easier to solve. 



VIII. CONCLUSION 



We have calculated the exact density of states and the typical localization length in a few interesting cases which 
are: the random electric potential model, the multiple disorder model, the random mass model in presence of a 
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constant magnetic field. Doing so, we used efficient algebraic methods to pinch the desired sectors of the Hamiltonian. 
For such quantities, the replica treatment is more straightforward thanks to the use of the replica index n viewed 
as a deformation parameter. The dynamical algebra of the effective Hamiltonian is u(l, 1) and allows to generate 
easily the groundstate sector from which stem the thermodynamical properties of the models. The vector space of 
this groundstate sector is identified as a supplementary representation of the dynamical algebra, labelled by n. An 
analog derivation is done in the supersymmetric formalism. 

We have also pointed out the existence of mappings between Hamiltonians describing different disordered models. 
In particular we mapped a random mass problem to random scalar potential problem by means of a non-unitary 
transformation. 

Those calculations can only be partially generalized to the calculation of the two-point Green function. The algebraic 
simplifications in the replica and supersymmetric cases we took advantage of in pointing out the dynamical algebra 
are still valid, as well as the use of the Laplace method. Yet the differential equations that come about are invariably 
the confluent equations of a six elementary singularities second order differential equation. Unfortunately this forbids 
exact calculations. Besides approximations can only be done safely in the case where there is a critical point and its 
simplifying properties, that is in the already well-know random mass case. 
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FIG. 1. Density of states and inverse localization length for the random mass model g x = 1. and m = 0. 
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FIG. 4. Density of states and inverse localization length for the random mass model in constant magnetic field g x = 1., 
m = 0.5 and 0=1. 
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TABLE I. Density of states and localization length of free massive Dirac particle. 
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TABLE II. Average integrated density of states and inverse localization length of a random mass Dirac particle in a constant 
magnetic field. M^(z) = J^(z) + N^(z), and J M and JV M are respectively the Bessel functions of first and second type. 
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